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Soft gluon resummation for heavy quark and dijet cross 
sections 
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Abstract: We discuss the resummation of threshold logarithms for heavy quark and dijet cross sec- 
tions in hadronic collisions. The resummed cross sections are presented at next-to-leading logarithmic 
accuracy in terms of anomalous dimension matrices which describe the factorization of soft gluons 
from the hard scattering. We apply our formalism to the calculation of the top quark production cross 
section at the Fermilab Tevatron. 



> 

o 
o 

ix. 



X 



1. Introduction 

Threshold resummations have recently been stud- 
ied for a variety of hadronic processes, in particu- 
lar heavy quark and jet production (for a review 
see Ref. Near threshold for the production 
of the final state in these processes one finds large 
logarithmic corrections which originate from soft 
gluon emission off the partons in the hard scat- 
tering. These logarithms can be resummed to 
all orders in perturbative QCD. The resumma- 
tion of the leading theshold logarithms in heavy 
quark and dijet production cross sections follows 
the earlier study of the Drell-Yan process H, ^ 
since these logarithms originate from the incom- 
ing partons and are thus universal. 

The resummation of next-to-leading loga- 
rithms (NLL) for QCD hard scattering and heavy 
quark production, in particular, was presented in 
Refs. Q 1^. At NLL accuracy one has to take 
into account the color exchange in the hard scat- 
tering. The resummation is formulated in terms 
of soft anomalous dimension matrices which de- 
scribe the factorization of soft gluons from the 
hard scattering. Applications to top and bottom 
quark production at a fixed center-of-mass scat- 
tering angle, — 90° (where the relevant anoma- 
lous dimension matrices are diagonal), were dis- 
cussed in Ref. Q. More recently the total cross 
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section was calculated for top quark production 
at the Tevatron |^]. Resummation for single- 
particle inclusive cross sections has been consid- 
ered in Ref. |]. 

The NLL resummation formalism for the 
hadroproduction of dijets has been presented in 
Refs. 1^, Jet production involves additional 
complications relative to heavy quark production 
because of the presence of the final-state jets. 



2. Resummation for heavy quark 
production 

We begin with the resummation formalism for 
heavy quark production. First, we write the fac- 
torized form of the cross section and identify sin- 
gular distributions in it near threshold. Then 
we refactorize the cross section into functions as- 
sociated with gluons coUinear to the incoming 
quarks, non-coUinear soft gluons, and the hard 
scattering. Resummation follows from the renor- 
malization properties of these functions. The re- 
summed cross section is written in terms of ex- 
ponentials of anomalous dimension matrices for 
each partonic subprocess involved. 

2.1 Factorized cross section 

We study the production of a pair of heavy quarks 
of momenta pi, p2, in collisions of hadrons ha and 
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hh with momenta Pa and pf,, 

ha{Pa) + hb{pb) -> Q{pi) + Q(p2) + X , (2.1) 

with total rapidity y and scattering angle 6 in 
the pair center-of-mass frame. The heavy quark 
production cross section can be written in a fac- 
torized form as a convolution of the perturba- 
tively calculable hard scattering H^j with par- 
ton distributions 4>fi/hi at factorization scale /x, 
for parton fi carrying a momentum fraction Xi 
of hadron h: 



dQ"^ dy d cos 6 



Ef dXa dXb . 2\ 



Q 



(2.2) 



where S = {pa + Pb)^, = {pi +P2Y, and 
we sum over the two main production partonic 
subprocesses, qq QQ and gg QQ. We note 
that the short-distance hard scattering H^j is 
a smooth function only away from the edges of 
partonic phase space. 

By replacing the incoming hadrons by par- 



tons in Eq. (2.2), we may write the infrared reg- 
ularized partonic scattering cross section, after 
integrating over the total rapidity of the heavy 
quark pair, as 



dQ"^ dcos( 



dz 



dxa dxb 



Xb 



(j)f/f{Xa,y^ 



X cr 



X (j>jr/jr{xb, e) S I Z 

Q 



ff-'QQ 



1 



z, ■ 



XbS ^ 
a,{^Ji^)\, (2.3) 



where the argument e represents the universal 
collinear singularities, and we have introduced 
a simplified hard-scattering function, ajj^QQ. 
The threshold for the partonic subprocess is given 
in terms of the variable z, 

z^^, (2.4) 
s 

with s — XaXbS the invariant mass squared of 
the incoming partons. At partonic threshold, 
Zma.x = 1, there is just enough partonic energy 
to produce the observed final state. We also de- 
fine a variable r = Zmin = Q^/S. In general, 
a includes "plus" distributions with respect to 



1 — z, which arise from incomplete cancellations 
near threshold between diagrams with real giuon 
emission and with virtual gluon corrections, with 
singularities at nth order in of the type 



n! 



ln"(l-z) 



TO < 2n — 1 , 



(2.5) 



If we take Mellin transforms (moments with 



respect to the variable r) of Eq. (2^), the convo- 
lution becomes a simple product of the moments 
of the parton distributions and the hard scatter- 
ing function a: 



dQ^ dcos9 



^~^f,f{N, ^l\e) cTjf^QQiN, Q/fi, e, a.ifi^)) (2.6) 



with the moments of the hard-scattering func- 
tion and the parton distributions defined, respec- 
tively, by 



4>iN) = 



dz z'^-^a{z) , 



dx 



<P{x) 



(2.7) 



We then factorize the initial-state collinear diver- 
gences into the parton distribution functions, ex- 
panded to the same order in Us as the partonic 
cross section, and thus obtain the perturbative 
expansion for the infrared-safe hard scattering 
function, a. 

We note that under moments divergent dis- 
tributions in 1 — z produce powers of In A^: 



ln"(l-z) 



-O (In 



m 

m — 1 



f 1 

N) 



■ ln"+i N 



The hard scattering function a is sensitive 
to soft-gluon dynamics through its 1 — z depen- 
dence (or the N dependence of its moments) . We 
may now refactorize (moments of) the cross sec- 
tion into iV- independent hard components H/l, 
which describe the truly short-distance hard scat- 
tering, center-of-mass distributions -0, associated 
with gluons collinear to the incoming partons, 
and a soft gluon function Sli associated with 
non-coUincar soft gluons. Note that the mass of 
the heavy quarks eliminates final state collinear 
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singularities in heavy quark production. This 
factorization is shown in Fig. 1 where // repre- 
sents a pair of hght quarks or gluons that pro- 
duce a heavy quark pair. Fig. 2 defines the 
soft gluon function, Sli, which describes the cou- 
pHng of soft gluons to the partons, represented 
by eikonal lines, in the hard scattering. / and 
L are color indices that describe the color struc- 
ture of the hard scattering. The hard-scattering 
function takes contributions from the amplitude 
and its complex conjugate, Hjl = h* l hi. 

Then we write the refactorizcd cross section 



as 



i 



cLtt' 

dQ^d cos 9 



IL 



iL 1 —,9,as{pL ) 



Q 



X Sli 1^-^,9, a,{^^) 
x^,/,(^iV,^,6) ^,-/,-(^7V,^,,^. (2.9) 

The center-of-mass distribution functions ip 
absorb the universal coUinear singularities of the 
initial-state partons in the refactorizcd cross sec- 
tion. They differ from the standard light-cone 
parton distributions by being defined at fixed 
energy rather than light-like momentum fraction. 
The moments of ip can then be written as prod- 
ucts of moments of and an infrared safe func- 
tion. 

2.2 Resummation 



Now, comparing Eqs. (2_^) and ( p.9| ), we see that 
the moments of the partonic heavy-quark pro- 
duction cross section are given by 



(N) 



^Pf/f{N,Q/^i,e) 
0///(iV,M2,e) 



IL 



IL I —,9,asi^I )j Sli ,6',Qs(m ) 



(2.10) 



where // denotes qq or gg, and we have used the 
relations ipq/q = '^q/q and (j)q/q = 4'q/q- Wc wiU 
now discuss the exponentiation of logarithms of 
A'' for each factor in the above equation. 



The first factor, (V'////0///)^ in Eq. (|2lo| ) 
is universal between electroweak and QCD hard 



processes, and was computed first for the Drell- 
Yan cross section The resummed expression 
for this factor is 



E'^f\N,Q-') 



X exp ■ 



dii' 



[iMsili")) 



-^ff{N,a,{^l'^))\] , (2.11) 

where 7/ is the anomalous dimension of the field 
of flavor /, which is independent of N , and 7// 
is the anomalous dimension of the color-diagonal 
splitting function for flavor /. i?(^)(as) is an 
A'^-independent function of the coupling, which 
can be normalized to unity at zeroth order. The 
exponent Eif) is given by 



dz- 



- 1 



,(1-2)'"S 



dX 

T 



[a.(Ag2)] 

+ B(f^ [as{{l-zr^Q^)] 

+ i^(^)[a,((l-z)2Q2)]| , (2.12) 

where the parameter ms and the resummed coef- 
ficients B'-^^ depend on the factorization scheme. 
In the DIS and MS factorization schemes we have 
ms = 1 and ms — 0, respectively. The function 
Aif) is given by 



(2. 



13) 



with Cf ^ Cf = (iVc - l)/(2A^c) for an incom- 
ing quark, and Cf — Ca — Nc for an incom- 
ing gluon, with Nc the number of colors, and 
K = (67/18 -7rV6)-5n//9, where Uf is the 
number of quark flavors. B^-'') is given for quarks 
in the DIS scheme by B^i^a,) = -SCpaJiAn), 
while it vanishes in the MS scheme for quarks 
and gluons. Finally, the lowest-order approxi- 
mation to the scheme-independent v'^-^^ is v^-^"^ = 
2CfajTr. 

Next, we discuss resummation for the soft 
function. The soft matrix Sli depends on TV 
through the ratio Q/ {Nfi); its TV-dependence can 
then be resummed by renormalization group anal- 
ysis. The product HjlSli of the soft fimction 
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Figure 1: Factorization for lieavy quark production near partonic tliresiiold. 





Figure 2: Tlie soft-gluon function Sli in wliicli tiie vertices cj, c*i^ link eikonal lines that represent the incoming 
and outgoing partons. 



and the hard factors needs no overall renormal- 
ization, because the UV divergences of Sli are 
balanced by construction by those of Hjl. Thus, 
we have [B 



H 



(0) 



IL 



AO) 



ic 



CD 



DL 



-L/ ~ {zI)lbSbaZs.ai, (2.14) 

where and S*'"^ denote the unrenormalized 
quantities, Zi is the renormalization constant of 
the ith incoming partonic field external to hi, 
and Zs,Li is a matrix of renormalization con- 
stants which describes the renormalization of the 
soft function, including mixing of color struc- 
tures, and which is defined to include the wave 
function renormalization for the heavy quarks. 



From Eq. (2.14), the soft function Sli satis- 
fies the renormalization group equation 



Sli — — (rg)LBS'B7 



'Sla(Ts) 



AI ■ 



(2.15) 



where is an anomalous dimension matrix that 
is calculated by explicit renormalization of the 
soft function. In a minimal subtraction renor- 
malization scheme and with e = 4 — n, where n is 
the number of space-time dimensions, the matrix 
of anomalous dimensions at one loop is given by 



g d 



(2.16) 



Using Eqs. (2.10), (2.11), and the solution 



of the renormalization group equation for the 
soft function, we can write the resummed heavy 
quark cross section in moment space as 



X TriiJ 



{2 [E^f^\N,Q^) 



7/.(a,(M'^))-7/./.(^,«.(/^''))] 



— ,e,asifi^] 
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i-Q/N 1 I 

X 5(l,0,a,(gViV2)) 
X Pexp / -^rs(a.(M")) , (2.17) 

^ J J 

where the symbols P and P refer to path ordering 
and the trace is taken in color space. 

We may simplify this result by choosing a 
color basis in which the anomalous dimension 
matrix Fg is diagonal, with eigenvalues A/ for 
each basis color tensor labelled by /. Then, we 
have 



s 



LI 



X exp 



dp, 

Q/N 



S 



LI 



1, 



iV2 



Xi{a,{f)) + \l{a,{f))] 



(2.18) 



Thus, in a diagonal basis, and with ^ = Q and 
R{f) normalized to unity, we can rewrite the re- 
summed cross section in a simplified form as 



xSLi{l,9,a,iQyN^)) exp 



E 



(2.19) 



where the exponent is 



4y)(7v,(?,g2) = - 



dz 



- 1 



1 - z 



+ g'^f\a,{{l~zr^Q^)] 

+ gi''^'f^\as{{l^zfQ%e]] , 

(2.20) 

with the functions gi, 172, and 53 defined by 



MLJf) 
93 



-Xi -Xl + li^'^^"> + . (2.21) 



2 2 

In the next section we present the soft anoma- 
lous dimension matrices for heavy quark produc- 
tion through light quark annihilation and gluon 
fusion; also we give NLO expansions of the re- 
summed cross section in both partonic produc- 
tion channels. 



3. Soft anomalous dimension matri- 
ces for heavy quark production 

3.1 Soft anomalous dimension for qq QQ 

First, we present the soft anomalous dimension 
matrix for heavy quark production through light 
quark annihilation, 

qiPa, To) + qijpb, n) Q{pi,ri) + Q{p2, r2) , 

(3.1) 

where the p^'s and r^'s denote momenta and col- 
ors of the partons in the process. 

We introduce the Mandelstam invariants 



s = {Pa +Pb) , 
h = {Pa ~PiT ~ iri^ , 
ui = {pb -PiY - "^^ 1 



(3.2) 



with m the heavy quark mass, which satisfy s -\- 
^1 + ui = at partonic threshold. 

For the determination of Tg we evaluate the 
one-loop graphs in Figs. 3 and 4. For the qq 
QQ partonic subprocess Fg is a 2 x 2 matrix, 
since the process can be described by a basis of 
two color tensors. Here we calculate the soft ma- 
trix in the s-channel singlet-octet basis: 

Cl — Csinglct ^rarb^r\r2 : 

C-2 = Coctct = {'^F)rbra{'^F)r2ri i (3.3) 

where the Tp are the generators of SU{3) in the 
fundamental representation. The result of our 
calculation in a general axial gauge is j|, ||, ||] 



Fii Fi2 
F21 F22 



(3.4) 



with 

Fii = 



a. 



Cf [Lp + ln(20^^) -I- ni] 



n \ti 



as Cf , f ui 



r22 — ■ 

X |cf 



4 In ( — ) - \n{2y/i^a'^b) - Lp -m 



.31n( ^ 1 - In 



+ Le 



(3.5) 
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Figure 3: Eikonal vertex corrections to Sli for partonic subprocesses in heavy quark or dijet production. 





Figure 4: Eikonal self-energy graplis for lieavy quark production. 



where is the velocity-dependent eikonal func- 
tion 



-/3 



1 - 2mVs A 1-/5 



In- 



1 + /3 



(3.6) 



with /3 = -^/l — Ani^/s. The gauge dependence 
of the incoming eikonal lines is given in terms 
of = {vi ■ ri)^/|np, with the gauge vec- 
tor, which cancels against corresponding terms 
in the parton distributions. Here the dimen- 
sionless velocity vectors are defined by = 
Qwf/v^ and obey vf = for the light incoming 
quarks and vf = 2m? /Q^ for the outgoing heavy 
quarks. The gauge dependence of the outgoing 
heavy quarks is cancelled by the inclusion of the 
self-energy diagrams in Fig. 4. 

Fs is diagonalizcd in this singlet-octet color 
basis at absolute threshold, /3 = 0. It is also diag- 



onalized for arbitrary [3 when the parton-parton 
cm. scattering angle \s = 90° (where ui = ii), 
with eigenvalues that may be read off from Eq. 



(3.5) 



We can expand the resummed heavy quark 
production cross section to any fixed order in per- 
turbation theory without having to diagonalize 
the soft anomalous dimension matrix. 

The NLO expansion for qq QQ in the MS 
scheme is 



~MS(1) 2 , \ B 

%Jqq(1 -z,m ,s, ti, Ml) = %-_QQ — 



TT 



X UCf 



"ln(l - z)' 


+ 


1 


1 - z 


+ 


1 - z_ 



Cf ( 81n( j-j -2-2ReL^ + 21n^-^ 
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Ca -3 In 



ReL, 



In 



(3.7) 



where the Born cross section. Similar 

qq^QQ 

resuhs have been obtained in the DIS scheme. 
These results are consistent with the approxi- 
mate one- loop results in Ref. . Analytical and 
numerical results have also been obtained for the 
expansions at NNLO [ij] in both pair inclusive 
and single-particle inclusive kinematics. 

3.2 Soft anomalous dimension for gg QQ 

Next, we present the soft anomalous dimension 
matrix for heavy quark production through giuon 
fusion, 

giPa, To) + g{pb, n) Q{pi,ri) + Q{p2, r2) , 

(3.8) 

with momenta and colors labelled by the pi's and 
Ti's, respectively. Here is a 3 x 3 matrix since 
the process can be described by a basis of three 
color tensors. The relevant graphs for the cal- 
culation of Ts are the same as in Figs. 3 and 4, 
where now the incoming eikonal lines represent 
gluons. 

We make the following choice for the color 
basis: 



C-Z — if {Tp)r2ri J 



(3.9) 



where d"'"^ and /"''^ are the totally symmetric 
and antisymmetric 5?7(3) invariant tensors, re- 
spectively. Then the anomalous dimension ma- 
trix in a general axial gauge is given by [Q, ^, || 



IT Z 



In {'^VaVb) 



In 



tlUl 



Lr — ni 



(3.11) 



We note that Ts is diagonalized in this basis at 
absolute threshold, /9 = 0, and also for arbitrary 
p when the parton-parton cm. scattering angle 
is 9 = 90° , with eigenvalues that may be read 



off from Eq. (|3.1l| ). 

Again, we may expand the resummed cross 
section for gg QQ at NLO and NNLO or 
higher. The NLO expansion in the MS scheme is 



„MS(1) /-, 2 , \ 

%Jqq(1 - z,m ,s,ti,ui) 



(J 



gg^QQ' 



ln(l - z) 



xUCa 



1 - z 



1 - z 



1 



{Cf - Ca)Rc Lp - [N^ - 1) In 



Nr. 



ft^uA 


-Cf 


\ m^s J 




--l)ln 


( i\u\ 


\vn?s ^ 




-ul) 



(3.12) 



where „ =, is the Born cross section, 

gg-*QQ ' 



B, 



QED 



1 



(3.13) 



ti Ui Am'^s 

Ui ti tiUi \ tiUi ^ 

and Kgg — {N^ — 1)^^ is a color average factor. 



4. Diagonalization and numerical re- 





"Fii ^ 




suits for qq tt 




T22 ^^31 
[ 31 4Ar^ i 31 i 22 J 


(3.10) 


In this section we give some numerical results for 
top quark production at the Tevatron. 

As we noted in the previous section, the soft 



where 



Fii 



[-Cf{L[, + 1) 



Ca 



r /' "1 

F31 = — In 1 — 



iln {AvaVb) + l-ni 



anomalous dimension matrices, F5, are in general 
not diagonal. They are only diagonal at specific 
kinematical regions, i.e. at absolute threshold, 
/3 = 0, and at a scattering angle 9 = 90°. In 
these cases the exponentiated cross section has a 
simpler form, Eq. ( 2.19| ), and numerical studies 



are easier to perform. Thus, the resummed cross 
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section at 6* = 90° for top quark production at 
the Fermilab Tevatron and bottom quark pro- 
duction at HERA-B were presented in Ref. 1^. 

In general, however, we must diagonalize Eg, 
i.e. we must find new color bases where Eg is di- 
agonal so that the resummed cross section can 



take the simpler form of Eq. (2.19). The diago- 
nal basis is given by C = CR where the columns 
of the matrix R are the eigenvectors of Eg. A de- 
tailed discussion of the diagonalization procedure 
has been given in Refs. ^. 

The resummed partonic cross section for the 
process qq QQ can be written as 



2 

E 

= 1/2 



dcos 9 



dSAfqq,ij{sA,0) 



(■5,54,6') 



(4.1) 



where sa = s + ti + m. The da''^^^^j{s, S4, 9)/dsi 
are color components of the differential of the 
Born cross section. The function fqq^ij is given 
at NLL by the exponential 

fqq,V = (^^P[Eqq,lj] = exp[Eqq + Eqq{X,, Xj)] , 

(4.2) 

where in the DIS scheme 



pDIS 
^qq = 



3CV 
2~ 



with llIq — 84/ {2m?) and A the QCD scale param- 
eter. The color-dependent contribution in the 
exponent is 



Eqq{Xi, Xj) 



+X' 



I * 



as 



A2 



A2 



(4.4) 



in both mass factorization schemes, where i,j — 
v'^i"> 12 (see Eq. (p^)) where 



1,2. Here A' = A 
the A's are the eigenvalues of the soft anomalous 
dimension matrix. The cutoff Scut in Eq. (4.1) 



regulates the divergence of as at low 54: Scut > 
S4,min = 2m'^K/Q. We choose a value of the cut- 
off consistent with the sum of the first few terms 



in the perturbative expansion |l^, in the 
range 30s4_min < Scut < 40s4^min, which corre- 
sponds to a soft gluon energy cutoff of the order 
of the decay width of the top, thus giving a natu- 
ral boundary of the nonperturbative region (see 
also the discussion in |l^). The central value 
in our cutoff range, Scut = 35s4,inin, corresponds 
to S4/(2to2) = 0.04 for m = 175 GeV/c^ and 
A = 0.2 GeV. 

After the diagonalization of the soft anoma- 
lous dimension matrix, the resummed partonic 
cross section is given at NLL accuracy by M 



a]^{s,m^) = - I dcos9 / 

J-l Js, 

1 da(°)(s,s4,e; 



s-2msi/^ 



dsA 



A1-A2P 



E^2 + |Ai-Enne^.^.ii 



dsA 



r^2 + |A2-rn|2 e^.^.- 



-^ELRe(e---) 

^ 2Rc ((Ai - En)(A2 - En)*e^«-i^)] .(4.5) 

To calculate the NLL resummed hadronic 
cross section we convolute the parton distribu- 
tions 4>i/hi for parton i in hadron h, with the 
partonic cross section: 



b ^1 

E 

q=u •' "^0 



(4.6) 



where a^^^^{TS,m ) is defined in Eq. (|4.5| ) and 
To = (m + + ScntY/S. 

In Fig. 5 we show numerical results for the 
tt production cross section at the Fermilab Teva- 
tron with VS = 1.8 TeV. We use the CTEQ 4D 
DIS parton densities ||l^. The NLO exact cross 
sections, including the factorization scale depen- 
dence, are shown as functions of the top quark 
mass. We also plot the NLO approximate cross 
section, i.e. the one-loop expansion of the re- 
summed cross section, calculated with Scut = 
and /i^ — rn^. We note the excellent agreement 
between the NLO exact and approximate cross 
sections. 
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Figure 5: The NLO exact and approximate and the NLL improved hadronic ti production cross sections in 
the qq channel and the DIS scheme are given as functions of top quark mass for pp colhsions at the Tevatron 
energy, \fS = 1.8 TeV. The NLO exact cross section is given for ^"^ = (soUd curve), 4m? (lower-dashed) 
and (upper-dashed). The NLO approximate cross section with Scut = is shown for /i^ — (lower 

dot-dashed). The NLL improved cross section, Eq. (4.7), is given for Scut = 35s4,miii (upper dot-dashed), 
30s4,min (upper-dotted) and 40s4,miii (lower-dotted). 



The scale dependence of the NLL resummed 
cross section is significantly reduced relative to 
that of the NLO cross section. In order to match 
our results to the exact NLO cross section we 
define the NLL improved cross section 

imp ICS NLO.approx , NLO, exact 

'^qq,\iad ~ '''gg,had ~ '^gg,had + ^qq.had ' 

(4.7) 

with the same cut applied to the NLO approxi- 
mate and the NLL resummed cross sections. 

We plot the hadronic improved cross section 
in Fig. 5 for [i? = m? along with the variation 
with Scut- The variation of the improved cross 
section with change of cutoff is small over the 
range 30s4,min < Scut < 40s4,niin- At m = 175 
GeV/c^ and = 1.8 TeV, the value of the im- 
proved cross section for qq tt with Scut/(2m^) = 
0.04 is 6.0 pb compared to a NLO cross section 
of 4.5 pb at /X = m. The contribution of the 
99 ^ QQ channel to top quark production at 
the Tevatron is relatively small. With its inclu- 
sion the total tt cross section at the Tevatron is 



7 pb. Gluon fusion is much more important for 
6-quark production at HERA-B Q. 

5. Resummation for dijet production 

Here we review the resummation formalism for 
the hadronic production of a pair of jets. Dijet 
production entails additional complications due 
to the presence of final state jets; otherwise, the 
formalism is analogous to that for heavy quark 
production. 

5.1 Factorized dijet cross section 

We consider dijet production in hadronic pro- 
cesses 

haiPa) + hipb) ^ Jl{Pl,Si) + J2(P2,<^2) + X , 

(5.1) 

where 6i and 62 are the cone angles of the jets. 
The partonic cross section is infrared safe since 
the initial-state coUinear singularities are factor- 
ized into universal parton distribution functions, 
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and the use of cones removes all the final-state 
collinear singularities. 

We study inclusive dijet cross sections at fixed 
rapidity interval 

^y^{l/2)\n[{p+p^)|{plp+)], (5.2) 

and total rapidity 

yjj = (l/2)ln[(p+ +p+)/(pr +P2 )] ■ (5.3) 

To construct the dijet cross sections, we define a 
large invariant, Mjj, which is held fixed. A nat- 
ural choice is the dijet invariant mass, Mjj = 
(pi +^2)^, analogous to for heavy quark pro- 
duction, but other choices are possible, such as 
the scalar product of the two jet momenta, Mjj — 
2pi ■ p2. We note that for large Mjj at fixed 
Ay we have a large momentum transfer in the 
partonic subprocess. The resummed cross sec- 
tion depends critically on the choice for Mjj. 
We shall see that the leading behavior of the re- 
summed cross section is the same as for heavy 
quark production when Mjj is the dijet invari- 
ant mass, while it is different for the other choice. 

The dijet cross section can be written in fac- 
torized form as 



da 



llahl,—^JlJ2 



\ - f dxa_ dXb 

■'•^ fajb=q,q,g 
2^ J, ,,2^ 



dMjjdyjjdAy 



X H 



Salt 



Mjj ^ Mj 



(5.4) 



where H is the hard scattering and the 0's are 
parton distribution functions. 

As for heavy quark production, we may re- 
place the incoming hadrons by partons in the 
above equation, and integrate over the total ra- 
pidity of the jet pair. Then, we can write the 
infrared regularized partonic cross section, which 
factorizes as the hadronic cross section, as 



dMjj dAy ~ J 



Af2 



AI 

z,—^,Ay,as{iJ,'^),Si,62] , (5.5) 



where the argument e represents the universal 
collinear singularities, and we have introduced a 
simplified hard-scattering function, (Tf, where f 
denotes the partonic subprocesses. The thresh- 
old region is given in terms of the variable z, 



(5.6) 



where, as before, s = XaX^S with S = [pa +pbY ■ 
Partonic threshold is at Zmax = 1 while the lower 

limit of Z is Zmin = T = MjjjS. 

By taking a Mellin transform of the rapidity- 



integrated partonic cross section, Eq. (5.5), we 
reduce the convolution to a simple product of 
moments. 



dr T^-i 



dMjj dAy 



f 

X af (iV, Mjjjyi, a,{p?), 8^,82) , (5.7) 



with (Tf(7V) and ^(A^) defined as in Eq. (|2.7|). 
We then factorize the initial-state collinear di- 
vergences into the light-cone distribution func- 
tions 4>f/fi expanded to the same order in as as 
the partonic cross section, and thus obtain the 
perturbative expansion for the infrared-safe hard 
scattering function, Ui. 

We can refactorize the cross section, as shown 
in Fig. 6, into hard components iJ/L, which de- 
scribe the truly short-distance hard-scattering, 
center-of-mass distributions ip, associated with 
gluons collinear to the incoming partons, a soft 
gluon function Sl/, associated with non-collinear 
soft gluons, and jet functions Ji, associated with 
gluons collinear to the outgoing jets. As for heavy 
quarks, here / and L denote color indices that 
describe the color structure of the hard scatter- 
ing. The refactorized cross section may then be 
written as 



dr T 



N~l ^P'/g/b 



f IL 



dMjj dAy 

Mjj 



X J, 



(/l) 



Mjj 

N, ,di 



,Ay S 



M 



■J.J 



Ay 



M 



J J 



J, 



if 2) 



N, 



Mjj 
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(5.8) 

This refactorization is similar to the heavy quark 
case except that we now include in addition out- 
going jet functions in order to absorb the final 
state collinear singularities. 

5.2 Resummed dijet cross section 



By comparing Eqs. (5/7) and (|5.8| ), we can write 
the refactorized expression for the Mellin trans- 
form of the hard-scattering function ai as 



af(iV) = 



IL 



X J( 



(/i) 



M 



J J 



Mjj 



Ay S_ 



f M, 



J J 



Ay 



Si] Ji 



if 2) 



(5.9) 



The results for the resummation of the uni- 
versal ratio "0/0 and the soft gluon function were 
given in the context of heavy quark production 
in Section 2. Explicit expressions for the soft 
anomalous dimension matrices for dijet produc- 
tion will be given in the next section. Thus, 
the only thing remaining to write down the re- 
summed dijet cross section are expressions for the 
resummation of the final state jets. 

The moments of the final-state jet with Mjj = 
{pi ^V-i? are given by ||] 



J, 



(/.) 



M 



J J 



,5,; =exp E[f^){N,Mjj) 



(5.10) 



with 

pMjj/N , , 

Elj^){N,Mjj)= / C(V^)(a,(M'^)), 

(5.11) 

where the first term in the series for C^y/)(as) 
may be read off from a one-loop calculation. The 
leading logarithmic behavior of the cross section 
in this case is not affected by the final state jets, 
so we always have an enhancement of the cross 
section at leading logarithm, as is the case for 
Drell-Yan and heavy quark cross sections. 

The moments of the final-state jet with M^ j = 



2pi ■ p2 are given by Eq. (5.1C) with 



E[f.^ {N,Mjj) 



1 2^^^ - 

dz 

1-2; 



+ B[^^ [asiil ~ z)M^jj)]] , (5.12) 
where the function A'-^'> is the same as in Eq. 



(2.13) and the lowest-order term in B^^^ may be 



read off from the one-loop jet function. The re- 
sults include a gauge dependence, which cancels 
against a corresponding dependence in the soft 
anomalous dimension matrix. The leading loga- 
rithms for final-state jets with M^j = 2pi ■p2 are 
negative and give a suppression to the cross sec- 
tion, in contrast to the initial-state leading-log 
contributions. 

Using Eqs. (U), ( pll ), and (|t^), together 
with the evolution of the soft function, we can 
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write the resummed dijet cross section in moment 
space as 

di{N)=Rlj^ expJ [e^^'\N,Mjj) 

I i—a,b 



-2 r" ^ [7/. {as{n'^)) - jf^f^ (TV, aAp"))] 

J fJ. 



exp<j E[f,)iN,Mjj, 
, i=i,2 



i-Mjj/N II 



use f-channel bases for the partonic processes in 
dijet production except for the processes qq — > gg 
and gg qq which are better described in terms 
of s-channel cofor structures. 

Since the full cross section is gauge indepen- 
dent, the gauge dependence in the product of the 
hard and soft functions, _ff must cancel 

the gauge dependence of the incoming jets, ip, 
and outgoing jets, J(f-)- Since the jets are inco- 
herent relative to the hard and soft functions, the 
gauge dependence of the anomalous dimension 
matrices Tg must be proportional to the iden- 
tity matrix. Then we can rewrite the anomalous 
dimension matrices as 

a. 



(f S )KL 



(5.13) 



E 

=a.h.l.2 



= {r^S')KL+OKL — 
IT 

C7(/,)i(-lni/, -ln2 + l-7rz), (6.4) 



with C/. — Cf {Ca) for a quark (gluon), and 

with Vi = {vi ■ n)^/|rip. Here the dimensionless 

This expression is analogous to Eq. ( |2.17D for and lightlike velocity vectors i)f are defined by 

= Mjjv^l^ and satisfy v} = 0. 

In the following subsections we will present 

(f) 

the explicit expressions for V gl for all the par- 
tonic subprocesses in dijet production, by eval- 
uating one- loop diagrams as in Fig. 3. The full 
anomalous dimension matrices can be retrieved 



heavy quark production except for the addition 
of the exponents for the final-state jets. 

We give explicit expressions for the soft anoma- 
lous dimension matrices, Ts, for the partonic 
subprocesses in dijet production in the next sec- 
tion. 

6. Soft anomalous dimension matri- 
ces for dijet production 

We consider partonic subprocesses 

fa {Pa,ra) + fb {Pb, u) ~> /l (pi, ^i) + /2 (p2, , 

(6.1) 

where the pi 's and 's denote momenta and col- 
ors of the partons in the process. To facilitate 
the presentation of the results for Ts we use the 
notation 

T = \n(^^^ +TTi,U = \n(^^^ +ni, (6.2) 
where 

S = [Pa+Pbf , t = {Pa-Plf ,U={pa - P2f , 

(6.3) 

are the usual Mandelstam invariants. We note 
that for the definition of the color basis we can 
choose any physical channel s, t, or u. We will 



from Eq. (|6.4| ) 



6.1 Soft anomalous dimension for qq qq 

First, we present the soft anomalous dimension 
matrix for the process 



qiPa,ra) + qipb,n) q{pi,ri) +q{p2,r2) 



in the t-channel singlet-octet color basis 
We find 0,117 1 



2CfT -^U 



(6.5) 



(6.6) 



(6.7) 



-2U -^(T~2C/). 

We note that the dependence on T is di- 
agonal in this i-channel color basis and in the 
forward region of the partonic scattering, T —^ 
—00, where Ts becomes diagonal, color singlet 
exchange is exponentially enhanced relatively to 
color octet. 
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6.2 Soft anomalous dimension for qq —i- qq 

and qq qq 

Next, we consider the process 

q{Pa,ra) +q{pb,rb) q{pi,ri) + q{p2,r2) , 

(6.8) 

in the t-channel singlet-octet color basis 

C2 = Sr^nSr^ri- (6-9) 

The anomalous dimension matrix is ^0 , |l7 
-^^{T + U) + 2CfU 2U 



2CfT 



(6.10) 



which also applies to the process 



q{Pi,ri) + q{p2,r2) -> q{Pa,ra) + q{pb,rb) . 

(6.11) 

Again, we note that the dependence on T 
is diagonal in this i-channel color basis, and the 
color singlet dominates in the forward region of 
the partonic scattering. 

6.3 Soft anomalous dimension for qq —> gg 

and gg qq 

Here, we present the soft anomalous dimension 
matrix for the process 

q{Pa,ra) + q{pb:n) -> .9 (pi , Ti ) + .g (pa , ^2 ) , 

(6.12) 

in the s-channel color basis 



Cl 
C2 
C3 



(6.13) 



We find 
I s' — — 

TT 






(6.14) 



U ~T 



2 ([/ - T) {U ~T) ^ (T + U) 



The same anomalous dimension describes also 
the time-reversed process [|[ |j 

gipi.Ti) +g{p2,r2) q{pa,ra) + q{pb,rb) ■ 

(6.15) 



6.4 Soft anomalous dimension for qg qg 

and qg qg 

Next, we consider the "Compton" process 

q{Pa,ra) +9{Pb,rb) q{pi,ri) +g{p2,r2) , 

(6.16) 

in the i-channel color basis 



-'rari'-'rt,r2 i 



Cl — Sr 

C2 = a-^^^{n)r,r. , 

C3 = tr^^%T^X^^^ . (6.17) 

The soft anomalous dimension matrix is pO| 

a. 



(6.18) 



X 



iCF + CA)T 







U 



CfT+^U 



2U 



U CFT+^^U 



2Na > 2 

which also applies to the process 

q{Pi,ri) +g{p2,r2) q{Pa,ra) +g{Pb,rb) ■ 

(6.19) 

We note that T appears only in the diagonal 
of the anomalous dimension matrix and that in 
the forward region (T — oo) the color singlet 
dominates. 

6.5 Soft anomalous dimension for gg gg 

Finally, we consider the much more complicated 
process 

9{Pa,ra) +g{pb,rb) -> g{pi,ri) -I- 3(^2,^2) • 

(6.20) 

A complete color basis for this process is given 
by the eight color structures |1(][ 



C2 
C3 



Taril jrbr2l , jr^ril frbr2l 



1 



C4 = Piira,rb;ri,r2) = g'^r^n'^r^rs , 

C5 = Pssira,rb;n,r2) = ^rf^-'-i-d'-'''--^ 

5 

C6 = PsAra,rb;r^,r2)^ln'''r^'^\ 

C7 = Pio+wi'>'a,rb;ri,r2) = ^(^r^r^'^rira 
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_A A \ _ _ fr^rxc fThr-iC 

C8 = P27{ra,n;ri,r2) = ^{Srarjrir2 

1 Q 

A A A A _ _.iraricirbr2C 



where the P's are t-channel projectors of irre- 
ducible representations of SU{3) ||l8|, and we use 
explicitly Nc = 3. 

The soft anomalous dimension matrix in this 
basis is M 



with 



3x3 = 



TT 



^3x3 03x5 
05x3 TsxS 



3T 
3U 
3{T + U) 



(6.22) 



and 

r5x5 = 

6T 


3U 
4 






TT 

3T 



(6.23) 



(6.24) 



-6C/ 



3U 
2 



3C/ 
2 



3C/ 
2 

6C/ 
5 



ST 



3f7 
2 











2t/ 
3 





-3U 


3[/ 






4 
5 



4[/ 



We note that the dependence on T is diagonal 
and that color singlet exchange dominates in the 
forward region of the partonic scattering, T 
—oo, where Ts becomes diagonal. This has been 
the case for all the processes we analyzed in t- 
channel bases; suppression increases with the di- 
mension of the exchanged color representation. 

The calculation of the eigenvalues and eigen- 
vectors of Fg/ for the gg ~* gg process is more dif- 
ficult than for the other partonic processes since 
we are now dealing with a 8 x 8 matrix; however, 
explicit results have been obtained [|o). The 
eigenvalues of the anomalous dimension matrix. 



are 










Ai = 


A4 


= 3 


TT 


T, 


A2 = 


As 


= 3 


as 

TT 





3 — {T + U), 

TT 



TT 



u 



(6.25) 



(6.21) The eigenvectors have the general form 



(3) 
0(5) 

0(3) 

.(5) 



z = l,2,3. 



i = 4. 



(6.26) 



where the superscripts refer to the dimension. 

(3) 

The three-dimensional vectors are given by 
i = 1,2,3. (6.27) 



= (3) 



Si3 



The five-dimensional vectors e^^\ Cg^"* and Cg'^^ 
are given by Q 



,(5) 



-15 " 




' " 


15 T 
2 U 




3 
2 


15 T 
2 U 







3 




3 3 T 

4 2U 


1 




1 



„(5) 
-6 



-15 





15 T 




2 U 


15 


15 T 


2 


2 U 



-3 
1 



(6.28) 



The expressions for e^^^ and eg^'' are long but can 
be given succinctly by 

"6i(AD 



(5) 



where 



b2i\d 
b3i\d 

1 



A, — — Ai 



(6.29) 



(6.30) 
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and where the biS are given by 



+ (40T^ - mu-^ - eor^jj + 52rt/^)A-] 

b2{X'i) = -|-[20r2-50[/T + 44t/2 
ZK 



+ (10r-5C/)A^] 



— — [40r^ - 64C/3 - 120r2f7 
2(7 /v 

130TC/2 + (20r2 + 13C/2 - 20TC/)Ay , 



h^{K) - -^(2T + 5C/-2A0, 



with 



K' = 2QT^ - 2WT + 21U^ 



(6.31) 



(6.32) 



7. Conclusion 

We have discussed soft gluon resummation at 
next-to-leading logarithmic accuracy for heavy 
quark and dijet production in hadronic collisions. 
We have constructed the resummed cross sec- 
tions in terms of exponentials of soft anomalous 
dimension matrices which describe the factoriza- 
tion of noncoUincar soft gluons from the hard 



"scattering! We have presented explicft resufts 



for the soft anomalous dimension matrices for 
the relevant partonic subprocesses. Numerical 
results have been presented for top quark pro- 
duction at the Fermilab Tevatron. Similar re- 
summations have been recently applied to direct 
photon [0, |, |l|l and -I- jet H, |^ production, 
as well as to calculations of transverse momen- 
tum and rapidity distributions JT2| , pT| for heavy 
quark production. 
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